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Abstract 

' The observed enhancement of pp-production near the threshold in radiative decays of J/ip and 

e^e~-annihilations can be explained with final state interactions among the produced NN system, 
where the enhancement is essentially determined by NN elastic scattering amplitudes. We propose to 
use an effective theory for interactions in a TV TV system near its threshold. The effective theory is similar 
to the well-known one for interactions in a NN system but with distinctions. It is interesting to note 
that in the effective theory some corrections to scattering amplitudes at tree-level can systematically 
be summed into a simple form. These corrections are from rescattering processes. With these corrected 
amplitudes we are able to describe the enhancement near the threshold in radiative decays of J/ ip and 
' e+e~-annihilations, and the pp elastic scattering near the threshold. 

CD 

r~| I It has been observed the enhancement of the pp-production near the threshold in various experiments. 

The enhancement has been observed near the threshold of the pp system in the decay J/ip — t- 'jpp by 
^SJ . BES[T] and in the decay of — t- K^pp and — )■ D^pp by Belle[2]. The enhancement also has appeared 
\ in the process e'^e~ — t- pp measured by Babar[3j. The enhancement of other baryonic system has been 
also seenjU [5]. In this work we focus on the enhancement of pp system. 

Because BES was the first to publish the result about the enhancement with rather high statistical 
I accuracy and precise information about the spectrum, many explanations for the observed enhancement at 
' BES exist. A class of explanations is that the enhancement is interpreted as the existence of a baryonium 
bound state[6] or a glueball below the threshold[7]. Another class of explanations is to take the effect 
of final state interactions into account. There are different ways to take final state interactions into 
account. One can use a complex 5- wave pp scattering length [8] or use a ii'- matrix formalism to include 
one pion exchange A more realistic way is by using potential models of NN interactions [9} [TO] . The 
observed enhancement in e"'"e~ — )• pp has also motivated theoretical studies [ini llSl [131 El ttSl IS]- It is 
I interesting to note that by taking final state interactions into account through potential models of NN 
interactions, the enhancement in J/ip — )• ^pp and e^e^ — )• pp can be explained simultaneously I15j . 
However, these models are in general complicated and contain several or more parameters which need to 
be fixed. It should be noted that with final state interactions the enhancement is predicted with the energy 
dependence of pp-elastic scattering amplitudes, where one employs Watson- Migdal approximation ^1 7j . 

A rather simple approach with final state interactions for explaining the enhancement has been give in 
[18j . The pp-elastic scattering amplitude there is determined through the rescattering mechanism similar 
to Watson-Migdal approximation, where one resums the multi-pion exchange or multi-rescattering of 
a NN system. Because the coupling of ttNN is well-known, the amplitude is completely fixed in this 
approach. It has been shown in [18j that the observed enhancement in the J/ip ^ ^pp and e+e^ — )• pp 
can be well explained. But, with the fixed amplitude one can not explain the pp elastic scattering near 
the threshold. In this work we make an attempt in the framework of effective field theories to give an 
unified explanation for the enhancement in the J/il^ decay and e~^e~ annihilation, and for the pp elastic 
scattering near the threshold. 
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An effective theory of NN interactions can be developed in analogy to the effective theory of NN 
interactions. The effective theory of A^A^ interactions has been proposed in [191 EOl El] and studied 
extensively in [20\ [2T\ I22j . With the effective theory the experimental data of A^A^ scattering of low 
partial waves near the threshold can be well described. In constructing such effective theories one makes 
an power expansion in the momentum near the threshold. The coefficients in the expansion characterize 
the properties of the A^A^- or A"A^ system, like scattering length oq, interaction range tq, etc.. 

There are distinct differences between the effective theory of A^A" interactions and that of A^A^ in- 
teractions. As for an effective field theory, one does not only need to construct the effective Lagrangian 
which gives tree-level amplitudes directly, but also one should be able to estimate relative importance 
of higher order corrections and to systematically calculate these corrections. In other word, one needs a 
power counting for the estimation of loop corrections. If the interaction of a system is characterized by a 

momentum scale A and the quantities like ao, tq, have the natural size, i.e., ao ~ A~^, a simple power 

counting like that for the well-known chiral perturbation theory can be used. However, for A^A^ systems 
it is not the case. It is well known from experiment that A^A^-systems have large scattering lengthes. 
This fact makes the simple power counting invalid. An important idea has been suggested in [21], where 
one can use the power divergence subtraction scheme instead of the minimal subtraction scheme which 
is commonly used, implemented to the effective theory. With the power divergence subtraction scheme a 
power counting can be established. This makes the effective theory well defined for A^A^ system with large 
scattering lengthes. In the case of A^A' systems, the scattering lengthes, according LEAR experiment[5] 
and model results [9], are around Ifm or smaller. They are much smaller than those of A^A^ systems. 
Therefore one can use the minimal subtraction scheme and hence the simple power counting for the 
effective theory of A^A" systems. 

Another difference is that a A^A^ system can be annihilated in to a multiple pion system and the pions 
can be real, while a A^A^ system can not be annihilated. The annihilation of a A^iV system into virtual- 
or real pions results in that the dispersive- and absorptive part of the A^A^ scattering amplitudes are of 
the same importance. In order to incorporate this fact some coupling constants in the effective theory 
of a A^A" system are complex numbers. In the effective theory of a A^A^ system the coupling constants 
are real. In this work we will first discuss the effective theory of a A^A" system and the A^A^ elastic 
scattering. Then we study the rescattering mechanism for the enhancement in the decay J /ip ^ 'jpp and 
the annihilation e~^e~ — )• pp. We will show that with our effective theory approach the enhancement in 
the decay J /ip ^ ^pp and the annihilation e~^e~ — )■ pp and the pp elastic scattering near the threshold 
can be well described. 

We consider the A^A^ scattering near the threshold: 

N{p,si) + N{-p,S2)^ N{k,s[) + N{-k,s'2), IpI = p = /3ym2+p2, (1) 

where p and k are three-momenta. The spins are denoted with s's. /3 is the velocity. Near the threshold, 
the momentum poi (3 approaches to zero. We are interested in the momentum region p < m^. An effective 
Lagrangian can be obtained by an expansion in p or /3. For this it is natural to use nonrelativistic fields 
to describe the nucleon A^. The nonrelativistic fields are given as 



■*= ' . x= . (2) 





The two-component field V'p(Xp) annihilates a proton(an anti-proton) and the field ipjiixp) creates a 
proton(an anti-proton). We denote the Pauli matrix acting in the SU{2) isospin space as r*(i = 1,2,3). 
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The TT fields are given as vr = rV* with vr^ = vr'^. At leading order the interacting part can be written as: 



with ~ 1.25 and -F-n- ~ 93MeV. cj with / = 0, 1 is the coupling constant in the channel, while dj 
with / = 0, 1 is the coupling constant in the ^Si channel. As discussed before, these coupling constants 
are in general complex because a NN can annihilate into vr's. One should keep in mind that the complex 
coupling constants here do not mean the violation of time-reversal symmetry. The complex coupling 
constants can be understood as the following: One can imagine that the effective theory is obtained from 
a perturbative matching of a more fundamental theory. In the more fundamental theory with the time- 
reversal symmetry amplitudes at tree-level are real, but they receive imaginary parts beyond tree-level 
because absorptive parts are nonzero at one- or more loop level. The imaginary parts of the coupling 
constants are from these absorptive parts in the matching. 

To clearly discuss the mentioned simple power counting for the above effective theory we first ignore 
the interactions of pion exchanges. Then in this case, the scattering amplitude at tree level is expanded 
in p, the leading orders are determined by the contact interactions given in Eq.(3) and are at 0{p^) if 
we take co,i and (io,i as constants. The tree- level contributions at higher orders of p starting at 0{p'^) 
are given by operators with derivatives in the effective theory. Therefore, at tree-level the amplitude 
is simply expanded in power of p and the power of p of each term is determined by the corresponding 
contact terms in the effective theory. However, this can be changed if we take loop-effects into account, 
i.e., the effects of scale-dependence of coupling constants. E.g., in the effective theory for a A^A^ system 
of large scattering lengthes the coupling constants corresponding to co,i and (io,i should be taken as at 
order of p~^ after including loop effects with power subtraction scheme, where one reasonably takes the 
renormalization scale as ~ p. In this way, a consistent power counting is established for a A^A^ system 
with large scattering lengthes [21]. 

As discuss before, it is expected that quantities characterizing interactions of a NN system have the 
nature size A. With this expectation the coefficients co,i and do,i in the effective theory scale like A~^ or 
(mA)~^ with m as the nucleon mass. They are at order of 0{p^) with the minimal subtraction scheme 
in dimensional regularization as shown in [21]. The scale-dependence of these coupling constants are 
suppressed by certain power of p. The loop contributions formed only by the contact interactions can be 
then estimated as the following: Each loop contributes a factor (mp/47r). Hence, the leading contribution 
is at p^ and comes from co,i and (io,i at tree-level. The next to-leading order is at p^ and comes from co,i 
and do,i at one-loop level. The contribution at p^ comes from co,i and do,i at two-loop level and from 
dimension-8 operators in 5C at tree-level, etc. In the above we discuss the simple power counting without 
interactions with pions. Since we consider the momentum region of p < m^, the interaction with vr is 
taken at the order of 0{p^). In Eq.(3) the leading interactions are given explicitly. At higher order of p 
operators with derivatives and operators for emission of more than one pion will appear. We notice that 
the power counting of loop-contributions through exchange of pions is slightly different because some 
loop contributions can be more strongly suppressed than the estimated by the simple power counting. 

With the interactions given in the above and also for our purpose, it is convenient to work with 
partial waves of the scattering. The scattering amplitude of Eq.(l), denoted as 7}(j>, k, si, S2, s'l, s'2) with 
the isospin I = 0, 1, can be decomposed with CG coefficients and harmonic oscillators into partial waves 
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Tlip,k,Si,S2,s'i,S2 



) = E ^4 iP/p)Yi'e', {k/k){l/2, si ■ 1/2, S2\s, S3) (1/2, ^ ; 1/2, 4 
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In the above the repeated indexes are summed and E is the total kinetic energy E = /m of the system. 
With the effective Lagrangian it is clear that the term with c/ will only contribute to T^oooo,/] (-^) with 
/ = 0, 1, respectively. 

With the effective theory it is straightforward to work out the scattering amplitude and various 
partial waves at tree-level to compare with experimental results of pp-scattering and enhancement near 
the threshold. In the comparison we perform a combined fit to fit the cross-section of pp-scattering and 
the enhancement in J /ip — t- 7p and e^e^ — t- pp, where the coupling constants co,i and do,i are taken as 
free parameters. However, as we will mention and show later, the theoretical results at tree-level can 
only fit the experimental data in a small region with E < 20MeV. Although the ^ jd.o.f . of the fit is 
close to 1, but the coupling constants can only be determined with the error from 40% to 100% or even 
more. In this work we will improve the situation by adding some corrections beyond tree-level. In the 
perturbation expansion of the effective theory some corrections appearing in n-loop level can be summed 
into a compact form. We will add these corrections to our tree-level results at tree-level and show that 
the improvement is significant. In the below we will study these corrections. This also illustrates some 
aspects of the effective theory. 





(6) 



Figure 1: The loop diagrams with the vertex of cq. The black dots denote vertices of the contact term. 



We focus at moment on the contribution only from cq by ignoring the contribution from 7r-exchange. 
It is straightforward to obtain the contributions from tree-, one-loop, etc. as given from the diagrams 
in Fig. la. The tree-level amplitude is just 47rco. The one-loop contribution is proportional to the loop 
integral regularized in (i-dimension |2 1] : 
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An interesting observation can be made for the above result. The integral is finite with d = 4. It has 
a pole at d = 3 corresponding to a power divergence of the integral. The power divergence subtraction 
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scheme is to subtract from the above contribution the pole contribution at d = 3. The subtraction 
introduces the renormahzation scale ^, hence a /U-dependence of the coupling constant. With this scheme 
one can show that for A^A^-interactions with large scattering lengthes a consistent power counting can 
be established [2 Ij by setting yu ~ p ~ m.,^. As discussed before, for A^iV-interactions scattering lengthes 
are much smaller than those from A^A^- interactions. Therefore we can employ the minimal subtraction 
scheme in which one subtracts the pole terms at d = 4. Because the above one-loop integral is finite 
at d = 4, no subtraction is needed. With d = 4 we have 47r/o = —imp = —im?/3. Inspecting the 
contribution from n-loop, one will find that the contribution is proportional to /q . In fact the sum of 
n-loop contributions forms a sum of a geometric series. The sum can easily be performed. We have then 
the exact the amplitude without vr-exchange as: 



'7(0000,0] (E) 



47rco 



Co 



(6) 



Expanding the above expression in co one can identify that the term with Cq"*"" comes from the n-loop 
diagram. Through the expansion one also sees that each loop brings a suppression factor p or /3 as 
indicated by the discussed simple power counting. The tree-level amplitudes obtained from the effective 
Lagrangian in Eq.(3) are at order of 0{p^). The corrections to them start at order of 0{p). By considering 
corrections from higher orders in p, the coupling constants in our effective theory will be generally depend 
on the renormahzation scale fi. With the power counting the /i-dependence is suppressed at least by 0{p). 






Figure 2: Fig. 2a represents the sum of rescattering amplitudes, where the nucleon- and antinucleon lines 
going through the cuts are on-shell. The narrow and long bubble represents the tree-level amplitude with 
the contact interactions and pion-exchange as indicated by Fig. 2b. 



Observing the above result one can realize that the loops in Fig. la calculated with the minimal 
subtraction scheme are the same in Fig. lb calculated only by taking the absorptive part of the loop 
integral lo, i.e., calculated by putting the A^A" in the loop on-shell indicated by the cuts in Fig. lb. Now 
we consider to include the contributions from vr-exchanges. The tree-level contribution can be represented 
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by the first diagram of Fig.2a where the bubble contains the vertex of cq and one-7r-exchange as indicated 
in Fig.2b. The contribution of Fig.2b is straightforward to obtain: 



7i£l)0,o](^) = 4vr 
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Co - -7:91 (2 - yLy] 



with the notations: 



9l 



( 9 A 



\2Fj ' ^ 2p2 

It is easy to find that the sum of bubble diagrams with cuts in Fig.2a can be performed. It is a geometric 
series. We have then the summed amplitude: 
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where • • • stand for other corrections, like the dispersive part of one-loop contributions formed with vr- 
exchanges. In the work we will neglect these corrections for our predictions. These corrections can be 
systematically studied. We will come back to the dispersive part in a future work, where a complete 
one-loop calculation will be performed. Physically the interpretation of Fig.2a and the amplitude in 
Eq.(9) is the following: The NN undergoes a multiple scattering process NN — t- NN — t- • • • — t- NN . 
Each scattering is due to the vertex with cq or exchange of one pion. Each pair of NN is on-shell in Fig. 2. 
We will call amplitudes for such a multiple scattering process as rescattering amplitudes. Similarly one 
can work out the tree-level and rescattering amplitude with 7 = 1: 
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Now we turn to the amplitudes with j = 1 and s = 1. In this case it is little complicated because 
the vrA^iV-interaction mixes amplitudes with difference i. The difference caused by exchanging one vr can 
only be ±2 . The summed or rescattering amplitude has to be expressed in a matrix form. We define 
the following matrix amplitude: 



M^j}{E) 



The tree-level results for M.[oj{E) reads: 



Tiiooi,i]{E), 7(1201,7] (^) 



(11) 



M'^l{E) = n 



V2gi (2 - 6y + {2y + 3y^)Ly 
V2gi (2 - 6y + {2y + 3y^)Ly) , g^[2 + 6y- (4y + 3y^)Ly] 



Mo + 2gl (2 - yL. 



(12) 



for 7 = 1 the matrix M^^{E) is obtained by replacing do with d\ and Qq with —gg/S. The summed or 
rescattering amplitude matrix from Fig.2 can be found as: 



M^ij{E)=Mf;i{E) 
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(47r)2 



(13) 



6 




Figure 3: The diagrammatic explanation of our approach. The dash line in Fig. 3a is the cut. 



It should be noted that the contact terms in the Lagrangian are involved only in the above discussed 
amplitudes. Later we will use these amplitudes and those of other partial waves at tree-level to compare 
with experimental results of the cross section of pp scattering near the threshold. 

Now we turn to the enhancement observed in J/ip ^ jpp and e~^e~ — )• pp. In [18| we assumed that 
in these processes a NN system is produced first, then the NN system undergoes a scattering before 
observed. This can be explained with Fig. 3. As illustrated in Fig. 3, The production of a pp system can 
be thought that at first step a nn system or a pp system is produced and then the NN system through a 
rescattering is converted into the observed pp system. We assume that the production of a NN system 
is a short distance process, i.e., the energy scale characterizing the production is much large than m-jr 
and the momentum near the threshold. Hence we can expand the vertex for the production in Fig. 3a in 
the momentum p or (3. At the leading order we can approximate the black vertex for the production in 
Fig.la with constant form factors. We denote the form factors as App^nn, respectively. At the leading 
order the NN- hence the final pp system can only be in the ^5*0 state. Then the decay amplitude can be 
written as: 



l + i 
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m?j3 
(4^ 



"^0000,1] 
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(14) 



The • • • stand for those partial wave amplitudes with £ > 1 and the corrections proportional to (3. With the 
approximation of the rescattering as discussed before, the pp amplitudes can be obtained from Eq.(9,10). 
We have then for the decay amplitude near the threshold as 



J/i/;— 



App + A-n 



(47r)2 '[0000,0]'^^^ 



1 -1 



+ 



App Ar, 



(0) . . 
(47r)2 ' [0000,1] l^^^ 



(15) 



Setting co,i = we recover our early results in [18] where in the corresponding amplitude the factor 
i should be replaced with —i because the absorptive part has been identified with a wrong sign. The 
physical results in |18j will be not changed with this replacement. In our approach the intermediate 
state can only be NN systems according to our effective theory. It is possible to add contributions from 
mesons as intermediate states in some phenomenological models as shown in [23j . 
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For the enhancement observed at Babar the form factors of proton are involved. These form factors 
are defined as: 



(p(Pi)p(P2)|J^|0) = n(pi) 



2m p 



v{p2), (16) 



It should be noted that near the threshold only the combination of the two form factors is involved. Near 
the threshold we have: 

{p{pi)PiP2)\m = -2mp + F2{q^)) ^^ar, + 0{(3^). (17) 



Similarly, we introduce the mechanism as in Fig. 3. We characterize the vertex of the A^A^ production 
at the first step with 
approach is given as: 



at the first step with constant form factors, denoted as G\/ and G\/ . Then the form factor with our 



Fiiq') + F,{q^) = + ^n(O) + ^ ^ii(l) (18) 
with Aii{I) as one matrix element of the matrix for the summed amplitudes as given in Eq.(13): 



A{I) 



(19) 



We will use the above formula for Fii^q^) + F2{q^) to describe Babar results. 

We will use our theoretical results to perform a fit by combining the BES data, Babar data in [3] and 
the data of the cross section of pp elastic scattering in pi]. For the BES data we use the measured results 
with BESS experiment [25]. In the fit we are able to determine the coupling constants in the effective 
theory, i.e., co,i and (io,i. These constants are complex as discussed before. We make two fits with 
different theoretical results in our approach. One fit is done with tree-level results, another is performed 
with NN scattering amplitudes in which part of partial wave amplitudes includes summed high-order 
corrections as indicated in Eq.(9-13) and in Eq.(15,18). The tree-level results for pp elastic scattering can 
be obtained from the effective theory in Eq.(3). In the fits we start to gradually add more data points to 
those nearest to the threshold until the /d.o.f larger than 1. We do not take these data points below 
the nn-threshold into account because we assume the isospin symmetry. Our fitting results are shown in 
Fig. 4 for BES and Babar and in Fig. 5 for the cross section. 

For the fit with tree-level results we are only able to describe the experimental data in the region with 
E <~ 20MeV. This can be seen from Fig. 4 and Fig. 5. Although the /d.o.f of the fit for experimental 
data in this region is around 1, but the coupling constants are determined with errors from 40% to 100% 
or even more. With the results containing the summed or rescattering partial waves the regions of data 
which can be fitted become larger as shown in Fig. 4 and Fig. 5. For BES data and Babar data the region 
is with E < 150MeV. For the cross section of the pp elastic scattering the region is with E < 70MeV. 
The results of the fit for the coupling constants in unit of GeV^^ are: 

Co = 123(15) - 22(ll)i, ci = -177(55) - 57(27)i, do = 32(31) + 17(60)z, di = 73(28) - 13(19)i. (20) 
The other constants appearing in J/ip ^ 'ypp and e'^e" — )• pp and x^/d-o.f. from the fit are: 



-\~ A-nn — 0.23(37), App A'fin — 2.22(72), 
Gfe) + gJJ) = 1.18(47), - = -0.11(20), 

xVd.o.f. = 50/63. (21) 
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™-2Mp(GeV) W-2Mp(GeV) 



Figure 4: The fitting results for the enhancement observed in J/-^ — )• ^pp observed by BES(left plot) and 
in e+e~ — )• pp by Babar (right plot). W is the total energy of the pp system. 




Figure 5: The fitting result for the cross-section of the elastic pp-scattering near the threshold 

In the above the coefficients App^nn are arbitrarily normalized. From the fitting results the coupling 
is not well determined. From Fig. 5 one can also see that without these contact terms in the effective 
theory the pp elastic scattering can not be described with pion exchanges. This is also true if we include 
rescattering effects in S'-waves. This is expected as discussed at beginning because without these contact 
terms in our effective theory the effects of a NN annihilation into virtual- and real pions are not included. 
These effects are important for NN scattering and should be not neglected. In our effective theory, the 
annihilation near threshold is described with these contact terms. 

From Fig. 5 we can see that in our approach for the pp elastic scattering near the threshold the 5- 
wave amplitudes are dominant. We can determine the scattering lengthes of S'-waves. In our notation 
the phase-shift and scattering length of a partial wave is given as: 

exp{2i6[ju,s,i]iE)) - 1 = *^'7[j«'s,/](-E^), pcot 5y«/s^/](^)) = h C'(/). (22) 

ovr C'ljWs,!] 
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Prom our fitting results we have the S'-wave scattering lengthes in unit of fm 



'^[0000,0] 
^[lOOl.O] 



1.80(19) + 0.32(17)i, a[oooo,i] = 2.61(81) + 0.84(41)i, 
0.48(46) - 0.25(89)i, a^iooi.i] = -1.07(41) + 0.19(29)i 



(23) 



These results are comparable with LEAR experiment [8j and model results^. The scattering lengthes 
are determined by the coupling constants in Eq.(3,9). We notice here that our numerical results of co,i 
in Eq.(20) are roughly twice larger than the corresponding coupling constants in the effective theory of 
a NN system [21]. But the determined scattering lengthes here are much smaller than the corresponding 
scattering lengthes of the A'^A'^ systems, because the relation between scattering lengthes and coupling 
constants is different in the two different effective theories. 

Before summarizing our study we briefly discuss our prediction relevant for the enhancement in 
J/ip ^ "YPp. The measured spectrum can be fitted with the decay amplitude as an 5- wave Breit-Wigner 
resonance form. The fitting done by BES gives the resonance mass around 1861MeV and the width 
r < 38MeV. The mass is near and below the threshold. We have observed that our formula in Eq.(15) 
can re-produce the shape of the S'-wave Breit-Wigner resonance form above the threshold only by tuning 
the parameters App^nn ■ If there is a resonance or structure below and near the threshold, one should 
also have roughly the same shape of the S'-wave Breit-Wigner resonance below the threshold. But, our 
amplitude in Eq.(15) below the threshold through analytical continuation has a cut near the threshold 
because of the L^-function in Eq.(8) appearing in the amplitude. Therefore, we can not conclude from 
our result that there exists a resonance or structure below the threshold, although our result agrees with 
experimental data above the threshold. 

To summarize: We have proposed an effective theory to study NN scattering near the threshold. 
Because scattering lengthes of AA^ system are rather small, we have implemented the standard minimal 
subtraction scheme to the effective theory to establish a power counting. The power counting is used 
to determine the relative importance of higher order corrections. We have found that certain higher 
order corrections represented as multiple rescattering can be simply summed. Using these rescattering 
amplitudes and assuming that the enhancement in J/ip ^ 'jpp and e^e~ — )• pp near the threshold of the 
pp system is due to final state interactions, we can simultaneously explain the experimental data of the 
enhancement and of the cross section of pp elastic scattering near the threshold. The S-wave scattering 
lengthes are determined which are comparable with existing results. Given the fact of the successful 
description of experimental data, the proposed effective theory needs to be studied in more detail. In 
the future we will study higher order corrections from loops formed through pion exchanges and the 
renormalization group of coupling constants in the effective theory. In this letter we have ignored the 
one-loop corrections to dispersive parts of scattering amplitudes, especially, the one-loop dispersive parts 
formed through pion exchanges. It has been found in AA-scattering that these parts can receive large 
corrections in [22]. It will be interesting to see if this also happens in A A scattering. This question can 
only be answered after our undergoing study of the complete one-loop correction has been done. 
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